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A general form for the surface roughness effects on the capacitance of a capacitor 
is proposed. We state that a capacitor with two uncoupled rough surfaces could be 
treated as two capacitors in series which have been divided from the mother capacitor 
by a slit. This is in contrast to the case where the two rough surfaces are coupled. 
When the rough surfaces are coupled, the type of coupling decides the modification of 
the capacitance in comparison to the uncoupled case. It is shown that if the coupling 
between the two surfaces of the capacitor is positive (negative), the capacitance is 
less (higher) than the case of two uncoupled rough plates. Also, we state that when 
the correlation length and the roughness exponent are small, the coupling effect is 
not negligible. 
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I. INTRODUCTION 


Nowadays widespread improvement in technological devices has offered operational ac¬ 
curacy in addition to miniaturization at submicron length scales. However, size effects due 
to miniaturization creates issues that are different from what experienced in macroscopic 
sized devices^. The fact of the matter is that by reducing the thickness of a thin film, the 
physical properties of the system vary. The principle issue here is that as devices tend to 
smaller dimensions, the coupling between their surfaces which is entangled to the existence 
of a cross correlation between the surfaces becomes more pronounced. Hence, by considering 
the fast growth in miniaturizing the devices, more attention should be paid on unignorable 
coupling effects. When the thickness decreases such that it gets comparable to the mean free 
path of the electron, the surface roughness comes in to play^. The surface roughness has 
been studied in the context of, e.g. electric conductivity^ - -, electron localization^, thermal 
conductivity^ - —, magnetization— - —, capacitance^, etching process^ 1 ^, leakage current^, 
wave scattering^ & shadowing effects^, surface growth^' 24 - & stochastic processes— - —, etc. 
A feature of surface roughness in the context of capacitance is linked to the fact that charges 
tend to accumulate on sharper areas. This statement proved adequate for Zhao et al. to 
show how the capacitance of a parallel plate capacitor increases when one of its plates gets 
rough22. Hence, it is obvious that the capacitance of a parallel plate capacitor which has 
two rough surfaces should be further modified. This issue provided the basis of this study; 
the encountered question here is wether a capacitor with two rough surfaces is equivalent to 
two capacitors in series each having one rough surface? It is instructive to state that in case 
of a non-coupled capacitor, if the capacitor is cut in half, the previous results obtained for 
a capacitor with one rough surface applies, see^. But when coupling exists, it is expected 
that the type of coupling between the two rough surfaces affects the physical properties of 
capacitors. In other words, one can expect to see different results obtained from a correla¬ 
tion or anti-correlation between surfaces. Note that when a surface is undergoing a growth 
process^ - —, the upper surface would not forget the previous information of the lower sur¬ 
face. Hence, the existence of correlation between the two surfaces is inevitable. Taking in 
to account the existing models in application to the configuration of the capacitor under 
consideration in this work, would lead to the elimination of coupling between the two rough 
plates of the capacitor. Hence, the coupling effects of the two rough surfaces is the center 


2 


of attention in this work. 


II. TWO BOUNDING COUPLED ROUGH SURFACES; LAPLACIAN 
SOLUTION 


Consider a parallel capacitor in which both surfaces are rough with a potential difference 
of V. The average distance between the two rough plates is d, where h\(x,y), h 2 (x,y) are 
the height fluctuations of the lower and upper plates respectively, see Fig. |T] The Laplace 
equation needs to be solved in order to obtain the electrostatic potential <E>(:r, y, z ) to provide 
basis for information about the physical properties of the system 

V 2 $(x,y,z) = 0, (1) 


where the boundary conditions for the potential obeys $(£, y,z — —d/2 + h\(x, y)) = 0 and 
<E>(x, y,z = +d/2 + h 2 (x, y)) = V. It is convenient to expand the boundary conditions using 
the Taylor expansion 


$(x,y,z = -d/2 + h 1 (x,y)) 
<F(x, y,z — +d/2 + h 2 (x, y)) 
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( 2 ) 


Assume that the roughness^ of the lower and upper surfaces are small compared to the 
average distance between the surfaces d, therefore it is instructive to utilize the perturbation 
expansion for the potential as 


<3>(x, y, z) = (x, y, z) + (x, y, z) + (x,y,z) + ■ ■ ■ . 


( 3 ) 


In this expansion, the nth order perturbed potential has an average as the order of 

(wi^/d) 71 . In order to find the perturbed potentials, we implement the techniques developed 
by Zhao et alM. These perturbed potentials individually satisfy the Laplace equation 

V 2 & n \x,y,z) = 0, (4) 


which is due to the fact that the terms have different orders of magnitude. It is worth stating 
here that the boundary conditions for the potential $ implies boundary conditions on each 
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Figure 1. Schematic of a parallel-plate capacitor with rough surfaces. The average distance between 
the plates is d, and the surface fluctuations are represented by hi(x,y),h 2 {x,y). 


of the perturbed terms which can be obtained by substituting Eq. (ED into Eq. ((2]) as 


M Uk pk 


k\ dz k 


M 


= ° , E it =C<5„,0, (5) 

t_n ,--d;2 tZo k a ~ ■- 

where M = 0,1,2, ■■■. It should be noted that for each value of M, Eq. (J5J) would give 
two recursive relations for the boundary conditions of perturbed potentials. By solving the 
Laplace equation for each perturbed potential, Eq. ([41) . with the consideration of its specific 
boundary conditions, Eq. (131) . the full solution for the potential ($) is obtained. Since the 
surface roughness is assumed to be small compared to the distance between the surfaces, only 
the first three terms of the perturbation expansion is considered and the rest are neglected. 
Hence, the zeroth order potential would be 

1/ / d\ 

( 6 ) 

and the first and second order potentials could be obtained making use of the Fourier integral 


$ ( 0 ) O ,y,z) = 7 ( z + 0 > 


as 


<5> {n) (x,y,z) = / df(q) 


,(n)/ .sinhgz (n) coshgz^ 


sinli 


+ A™(q) 


cosh y 


exp(—iq L .p)d 2 q, 


where we have n — 1,2, p = (x,y) and 


^(q) = “^(^(q) T^i(q)), 

^± 2) (q) = ^ f ^2(q-q')T/ii(q-q , )](^2(q , )-^i(q , ))9 , coth^dV 

+ ~ki I [^(q-q')±^i(q-q')](^2(q / )+ ^i(q , ))9 , tanh^d 2 g, 


(7) 


( 8 ) 
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Note that h( q) = jyjt f h(p)e zclp d 2 p is the Fourier transform of h(p). 

According to Eq. (I6j) . it could readily be noticed that the zeroth order potential is the 
indication of a smooth parallel-plate capacitor. For the first order potential, as shown in 
Eq. m an d Eq. (J8]), the term (A^ 1 ) contains information about the height fluctuation hi , 
where i = 1,2—. For the second order potential, as indicated by Eq. (J7| and Eq. (JHD , the 
term (A±^) contains information about the product of the height fluctuations hihj, where 
we have i, j = 1, 2. Strictly speaking, the nth order potential consists of terms including the 
product of the n height fluctuations ■ ■ ■ hi n , where we have ■ ■ ■ ,i n = 1, 2. 


III. EXPRESSIONS FOR THE ELECTRIC FIELD AND CAPACITANCE 


we proceed to find the electric held of our configuration. Since the height fluctuations 
of the surfaces are considered as a random held obtained by a distribution function, e.g. 
Gaussian, the potential <f> would accordingly be a random quantity, where its average at 
every point inside the capacitor is of our interest. By knowing the potential, the electric 
held could be obtained. Hence, for the electric held we have 


E(x, y, z ) = -V$ = -V$ (0) - V$ (1) - V$ (2) + • • • 

2 
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+ A H (q) ^ ) exp(-iq ,p)d 2 q. 
cosh y 


(9) 


Note that the electric held is a random vector held, and e 3 is the unit vector in the £ 
direction. The average of the electric held is also an interest in this work. 

In order to obtain the total charge on a rough-surface of a capacitor we recall the fact that 
the amount of charge accumulated on a typical point of the surface is proportional to the 
electric held normal to that point. Thus, the total charge accumulated on the lower {i = 1) 
and upper (i = 2) rough surfaces is given by Qi = f E,.^ ds, with surface normal vector 
nj = (—1)®—^=h=|L=. In order to obtain the charge of the plates ( Qi), information about 
the plates electric held (E* = E(x, y,z = (—l) l d/2 + hi)) is essential. Note that in order to 
keep the calculations a bit simple, since the roughness is smaller compared to d , one can 
write Ej ss E (x,y,z = ( — l) l d/2). Hence, the total charge accumulated in the capacitor is 
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readily obtained. Interestingly the accumulated charge is non-zero, Qi + Q 2 ^ 0. This is in 
contrast to Gauss’s law for the total charge accumulated on the surfaces of a capacitor; as it 
speaks for itself it should be zero. Hence, something is odd here. This inconsistency is born 
from the two different kind of assumed approximations: (%) only considering the first three 
terms of the perturbation expansion in Eq. (jUJ), and (ii) using the approximated version 
for the plates electric held, E f « E (x,y,z = (—l)*d/2). This discrepancy could easily be 
overcome by taking Q = (Q 2 + |Qi|)/2, as magnitude of the charge accumulated on each 
surface. Hence, the random capacitance is readily obtained 

C = Q/V. (10) 

It is worth noting here that the height fluctuations of the plates behaves in a random 
manner, so the ensemble average is of interest. Hence, we take the average functional form 
of the capacitance together with the electric held and potential. The fact of the matter 
is that in the process of obtaining the ensemble average of $, E and C, the height-height 
correlation function of the surfaces comes in to play. In the case of a capacitor with only 
one rough surface it was shown by Zhao et al. that for obtaining the ensemble average, the 
autocorrelation of a rough surface appears only in the second order perturbed term see 
re£2£. This motivates the study of the cross correlation effects for the case of a capacitor 
with two coupled rough surfaces. We intend to show that in the averaging process, the cross 
correlation effects between the two rough surfaces also shows itself in <fh 2 h 
There are two basic statistical properties of importance for describing random processes 
(or random fields). The correlation functions Rij{ ri,r 2 ) = (hi(ri)hj(r 2 )), and the spectral 
density functions SV,(q) = (h t ( q) hj (— q) ), where (■ ■ ■) denotes the ensemble average 

over possible roughness configurations. The parameter A denotes the area of the projected 
plate of the capacitor on the x-y plane, with i = 1, 322r— . 

In the present work, the joint distribution function of the height fluctuation is considered 
Gaussian with homogeneous and isotropic surfaces. It is well known that for a homogeneous 
and isotropic rough surface the correlation functions Rij{ r 1; r 2 ) = Rij {liq — r 2 |) and the 
spectral density functions Sij( q) = Sij(q) are real functions where q — |q| — . For a homo¬ 
geneous and isotropic rough surface due to the Wiener-Khintchine theorem^ -35 the spectral 
density is the Fourier transform of the correlation function. In addition, one can show 

{h i (q 1 )h j (q 2 )) = ^p<5(qi + q^K^qOM-qi))- ( n ) 


6 




To comply with the aims of this work it is essential to obtain the three main parameters 
of a capacitor; the average ensemble of the potential, electric held, and capacitance. Making 
use of Eqs. dSD-anD, and keeping in mind that due to the Wick’s theorem for a Gaussian 
distribution, the ensemble average of the product of any odd number of h( q) is zero, we can 
obtain 

($) « <f> (0) + (<f> (2) ) 

= 7 (^ + + ^ 2 JAd / (^ 2 ( q )| 2 ) “ (Mq)| 2 »9 cothqdtfq, (12) 

for the potential, and 

<E} = -(V4)«-^e 3 -(2x) 2 )Te3, (13) 

for the electric held, and 

= i + vg- p + -^r / ((Mq^ + t^q)! 2 })? 2 ^}, (14) 

for the capacitance. Where 

V = f (dM q )| 2 ) + (Mq)| 2 ) - (hi(q)h 2 (—q)) - (h 2 (q)hi(-q ))^q coth qdd 2 q. (15) 

It could readily be seen in Eq. (fT2l) that the autocorrelation and coupling effects show 
themselves in the term for (dh 2 '). In addition, in the second term of Eq. (fT2lh the contri¬ 
bution of the coupling and autocorrelation are the same in the sense of order. So when the 
contribution of coupling is as the same order of the autocorrelation, the effects of coupling 
is not negligible. Note that the last term on the RHS of Eq. (1T21) is just a constant term. 

In order to obtain the ensemble average of the electric held in Eq. (fT3ll we use the fact 
that (• • •) and V can commute with each other. The ensemble average of the electric 
held has no component on the x-y plane, this is expected for homogeneous and isotropic 
surfaces. Moreover, (E) is a uniform electric held as in the case for the parallel-plate 
capacitor. For the capacitance in Eq. (fT4l) the coupling effect enters in the same order as 
the autocorrelations of the surfaces which is similar to that for the potential and electric held. 
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IV. EFFECTS OF SURFACE COUPLING ON ELECTRICAL 
PARAMETERS 


To discuss our results three case studies are carried out; in case of a single rough sur¬ 
face, assuming that the upper surface is rough, the only non-zero term in Eqs. (fT2l) - (fT4l) is 
(|h 2 (q)| 2 ). This conclusion resembles the results obtained by Zhao et alM where the capaci¬ 
tance increases compared to a parallel-plate capacitor due to accumulation on sharp places, 
see Fig. |T| For the case where two uncoupled rough surfaces exist, the two coupling terms 
(h 1 (q)/r 2 (—q)) and (h 2 (q)hi(—q)) in Eqs. (fT^ - dTH) disappear. Thus, the presence of two 
uncoupled rough surfaces would cause an increase in the ($) , (E) , and ( C) compared to a 
capacitor with only one rough surface. For the case where rough surfaces are coupled, both 
autocorrelation and coupling terms count. In this case, depending on the functional form 
of coupling between the two surfaces we experience an increase or decrease in the values of 
( < F), (E) , and ( C ) in comparison to two uncoupled rough surfaces. 

The discussions carried out before this point are general, in a sense that any statement born 
out in application to the physical quantities of the capacitor up to now disregards the type of 
roughness of the surfaces. But, since most of the rough surfaces in nature may be considered 
self-affine, it is best to treat them accordingly, see — and references therein. An analytic 
model for the roughness power spectrum which has the proper asymptotic limits that allows 
calculation of roughness effects is given in —. This model has the Lorentzian form 

2v a w 2 e 

} (2vr) 3 (l + aq 2 e) W 

which proved to be consistent with observed data^. The roughness exponent a indicates 
the degree of roughness irregularity at short length scales r < £, where £ is the correlation 
length —. The parameter a in Eq. (fT6l) is introduced in a piece-wise form as 




| h\ l ~ ( 1 + a <lU 2 ) 1 0 < a < 1 

17jln(l + aq 2 £ 2 ) a = 0. 


(17) 


Note that the parameter q c = vr/ao is the upper limit for the frequency in the Fourier space, 
where ao is of the atomic order—. It could be noticed from Eq. (1T6l) that for a self-affine 
surface, the three parameters a, £, and w, specify all information about the surface. However 
these parameters could be different for various surfaces. 



In order to show the effect of coupling between two rough plates we suppose the surfaces 
self-affine. Consider two uncoupled self-affine rough surfaces characterised by ( 07 , £ 1 , uq) 
and (q? 2 , £ 2 j W 2 ) • The spectral density of each surface is obtained by Eq. (fTUl) where the 
coupling terms disappear as discussed earlier. If we take the roughness of the two surfaces 
as Wi = w 2 = w, and the roughness exponent of the two surfaces as aq = a 2 = a where 
w and a are constant, the only variable that remains is the correlation lengths, £i i2 . In 
this stage it is convenient to perform a change of variable and take A = w/d, L\ >2 = ^ 1 , 2 /d 
and q' = qd — . By substituting the spectral density functions of each self-affine surface in 
Eq. (US]), the electric held is obtained 


(E h 


2 ft.1 

l + A/ 

i=i J « 


q' 2 coth q' 


(1 + q' 2 L 2 / 2a) 1+c 


dq 


(18) 


where the index X stands for two independent or uncoupled surfaces, and Eq = V/d is the 
electric held of a parallel-plate capacitor. Similarly the capacitance is 


(Ch 


(E)x , A 2 




<?' 


,/3 


i =1 


/ 0 (1 + q ,2 L 2 /2a) 1+a 


dq 


(19) 


Co Eq 

where Cq = eo A/d is the parallel-plate capacitance. Figure [2] demonstrates the dependence 
of (C)x/C'o on the normalized correlation lengths Li j2 . Figure[ 2 ]is plotted for three diherent 
values of a with A = 0.01. It could be deduced from Fig. [2] that for a hxed normalized 
correlation length Li^, as the roughness exponent a decreases, the ratio (C)x/Co increases. 
In addition, for a hxed roughness exponent a, as Li or L 2 increases, the ratio (C)i/Cq 
decreases to unity. 

In the case of two coupled rough surfaces, the discrepancy lies in the fact that there is a 
non-zero cross-spectral density. In general, cross-spectral density (/ii(q)/i 2 (—q)) is equal to 


Ti2<q)^ (IMq)l 2 XMq)l 2 ). 


( 20 ) 


where 712(q) called the coherence function is a complex function located in the unit circle of 
the complex plane^ 1 ^. For two homogeneous and isotropic surfaces, the coherence function 
(712 (q) = 721 (q) = 7(9)) is real belonging to the domain [—1,1], Here, we suppose the simple 
case 7 (q) = —1 which indicates the negative cross-correlation between surfaces. Hence 
implementing the same change of variables for the roughness and correlation lengths, we 
notice that ((E) c - (E) x ) /E 0 = ((C) c - (C) x ) /C 0 equals 

2 A 2 q ' 2 coth q' dq' 

[(1 + q' 2 L 2 / 2 a)(l + q' 2 L 2 / 2 a )}^+ a )/ 2 ’ ( ) 
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Figure 2. Three Surface plots showing the ratios of the capacitances ( C)x/C ® for the case of two 
uncoupled rough surfaces in terms of the normalised correlation lengths L 12 (= £ 1 , 2 /d) . From top 
to bottom the roughness exponent has been taken equal to 0.3,0.5, 0.7, with A = w/d = 0.01. 

where the index C indicates two coupled surfaces. Figure [3] shows the coupling effect for 
both the electric field and capacitance Eq. CD as a function of the normalized correlation 
lengths L\ 2 for various roughness exponents a at fixed A = 0.01. Also, it could be deduced 
that for a fixed normalized correlation length L lj2 , as the roughness exponent (a) decreases, 
the capacitance of the coupled capacitor increases significantly in comparison to a smooth 
parallel-plate capacitor. Moreover, for a fixed a, as L\ or L 2 increases, a decrease is seen 
in (E)c and (C)c of the coupled capacitor. In other words, as the normalized correlation 
lengths prolongs, the system would tend to the case of an uncoupled rough capacitor. It 
could readily be noticed from Fig. [3] that the coupling effect is more efficient when the 
normalised correlation length and the roughness exponent are small. 

V. CONCLUSIONS 

In this work the Laplace equation was solved for a parallel-plate capacitor with two rough 
surfaces. Assuming that the roughness of both surfaces are small compared to the average 
distance between the surfaces, the perturbation of the electric potential was substituted in 
the Laplace equation. Since the fluctuations of the surface height was considered very small 
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Figure 3. Three surface plots illustrating deviation of the coupled capacitance from the uncoupled 
capacitance ((C)c — (C)x)/C' 0 ' in terms of the normalised correlation lengths L\^{= £ 1 , 2 /d). 
From top to bottom the values of the roughness exponent is taken equal to 0.3,0.5, 0.7, and we 
have taken A = w/d = 0.01. The same goes for the deviation of the corresponding electric fields 
(( E)c-(E) x )/EW . 


compared to the width of the capacitor, only up to the second order terms were kept. 
Solutions to the Laplace equation giving the electric potential of a capacitor with two rough 
surfaces (either coupled or uncoupled) leads to the conclusion that the electric held and 
consequently the capacitance increases compared to the case of the smooth parallel-plate 
capacitor, see also^H. This could be explained by the fact that electric charges accumulate 
on sharp places of a substrate where in the particular case studied here is the summit and 
foothill of the height fluctuations on the surface. The comparison of a capacitor with two 
uncoupled rough surfaces with another capacitor with two coupled rough surfaces showed 
that the increase and decrease of the capacitance for the two cases depend on the sign of their 
cross correlation. If the two rough surfaces are correlated/anti correlated, the capacitance 
is decreased/increased (Eq. ([IT]) ) in comparison to the capacitor with two uncoupled rough 
plates. Note that when the normalised correlation length and the roughness exponent are 
small, the coupling effect is not negligible, see Fig. [3] 

The model considered in this work was based on the fact that in a rough surface capacitor, 
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the coupling between the surfaces causes a deviation in its capacitance in comparison to 
what obtained by entering a slit in between its plates. In other words, the slit creates 
two capacitors in series (each with one rough surface) where their equivalent capacitance 
varies from the initial capacitance that had two coupled rough surfaces. The reason for this 
most interesting deviation is linked to the coupling between the two rough surfaces. To 
be more precise, the deviation would be most pronounced when the coupling between the 
rough surfaces is strong. The term strong comes from the fact that the cross correlation of 
the two rough surfaces is as of the same order of the height-height auto correlation of each 
rough surface. The contribution of a strong correlation towards the physical parameters of a 
capacitor could rise up to twenty percent or even more. Hence, considering a capacitor with 
two coupled rough surfaces as two capacitors in series where each has one rough surface may 
not be the best assumption. Hence, we understand now that the coupling effects should be 
taken more seriously into account as devices tend to miniaturize down to submicron ranges. 
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